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1.  Introduction 


This  thesis  examines  the  performance  of  coherent  phase-shift-keyed  (PSK) 
direct-sequence  spread-spectrum  (DS/SS)  communications  in  a  specular  mul¬ 
tipath  fading  environment.  In  a  previous  discussion  of  this  particular  topic 
[1],  the  authors  derived  an  expression  for  the  average  probability  of  error  of 
a  correlation  receiver,  based  on  the  characteristic-function  of  the  multipath 
interference,  given  an  arbitrary  number  of  paths  with  both  deterministic  or 
random  gain  coefficients.  The  channel  model  considered  here  is  a  direct 
extension  of  that  presented  in  [1]  with  some  important  constraints  and  mod¬ 
ifications.  First,  there  always  exists  a  nonfaded  component  in  the  zeroth 
path  signal,  i.e.,  the  direct  path  plus  multipath  interference  with  relatively 
small  delay  in  comparison  to  the  direct  path.  Second,  the  range  of  the  path 
delays  is  restricted  such  that  intersymbol  interference  is  due  to  the  same  pair 
of  adjacent  data  bits  for  each  of  the  reflected  paths  (case(i.)  in  [1]).  Third, 
for  the  sake  of  brevity,  only  deterministic  gain  coefficients  are  considered. 
Finally,  the  assumption  that  the  path  delay  must  be  at  least  one  chip  period 
is  removed  and  consideration  is  given  to  delays  less  than  a  chip. 

As  in  [1],  the  average  probability  of  error  is  approximated  using  the 
characteristic-function  method.  Specifically,  the  characteristic  function  of 
the  multipath  interference  is  numerically  integrated  using  Simpson's  method. 
This  technique  is  extended  also  to  compute  the  probability  mass  function  of 
the  receiver  output. 
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The  thesis  is  presented  in  the  following  fashion.  Section  2  discusses  the 
system  and  channel  models.  Section  3  provides  a  derivational  analysis  of 
the  average  probability  of  error  for  any  desired  range  of  path  delays  within 
model  constraints.  In  Section  4,  an  analysis  is  provided  for  the  performance 
of  channel  coding  applications.  Then,  the  models  are  numerically  evaluated 
for  the  PSK  DS/SS  system  in  Section  5.  Finally,  the  conclusions  of  these 
results  are  presented  in  Section  6.  Following  the  text  are  two  appendices. 
Appendix  A  links  the  characteristic-function  method  to  an  expression  for 
the  average  probability  of  error.  Appendix  B  contains  the  computer  code 
used  to  produce  the  numerical  results. 
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2.  System  and  Channel  Models 


Transmitter 

As  presented  in  [1],  a  channel  user  starts  with  a  data  signal  b(t).  a  se¬ 
quence  of  rectangular  pulses,  each  of  duration  T  .  where  the  /- th  pulse  has 
amplitude  bi  for  IT  <  t  <  (I  4-  1)J.  The  data  ( bi )  is  viewed  as  a  sequence  of 
mutually  independent  random  variables,  each  with  equal  probability  of  being 
-hi  or  -1.  The  signal  is  then  spread  by  the  waveform  x|/).  given  by 

x(f)  =  Xiv(t  —  iTc).  iT:  <  t  <  (i  -h  1  )TC  (1) 

where  x,-  is  the  binary  signature  sequence  and  c(s)  is  a  rectangular  chip 
waveform  satisfying  v(s)  =  1  for  0  <  $  <  T:  and  u(s)  =  0  otherwise  [2]. 
T:  represents  the  chip  duration  for  the  transmitted  signal.  The  processing 
gain  A"  is  defined  as  the  ratio  between  the  period  of  the  data  signal  and  the 
chipping  period:  i.e..  -V  =  T/Tc.  Thus,  the  signature  sequence  (x,)  has  a 
period  of  .V  also.  Considering  the  transmitter  power  P  and  the  phase  angle 
of  the  PSK  modulator  0.  the  created  spread-spectrum  signal  u{t)  can  be 
expressed  as 

u(t)  =  \/2Pb(t)x(t)  exp(jO).  (2) 

Thus,  with  carrier  frequency  f:  .  the  transmitted  signal  for  a  PSK  DS/SS 
system  is  given  by 

s{t)  =  Re[u(t)  exp(j'2~  f:t)].  (3) 
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Channel 

The  output  signal  y(t)  produced  by  the  specular  multipath  channel  con¬ 
sists  of  one  direct  path  and  L  delayed  paths  and  may  be  modeled  as  the  sum 
of  attenuated,  phase-shifted,  and  delayed  versions  of  the  input  signal  s(t). 
Thus,  as  in  [1],  it  may  be  expressed  as 

l 

y(t)  =  Rei9nu{t  ~  rri)exv[j2-fc(t  -  t»)]}  (4) 

n— 0 

where  the  complex  gain  coefficients  gn  are  given  as 

go  =  3  +  -7o.4o exp(j0,)  and  gn  =  -yn.4nexp(;0n).  1  <  n  <  L  (5) 

for  the  nth  delayed  path.  Here.  */n.4„  and  9n  represent  the  attenuation  and 
phase  shift  incurred  from  the  reflection  of  the  delayed  paths,  respectively, 
while  rn  is  the  nth  path  delay  relative  to  the  zeroth,  or  direct  path,  assuming 
that  t0  =  0.  The  parameter  3  signifies  the  gain  coefficient  for  the  nonfaded 
component  of  the  direct  path.  Throughout  the  analysis,  it  is  assumed  that 
3  =  1  and  the  random  variables  .4n  are  normalized  to  have  second-order 
moments  equal  to  1.  Also,  may  be  interpreted  as  the  relative  power  of 
the  nth  delayed  path  as  compared  to  the  nonfaded  component  of  the  direct 
path. 

As  previously  mentioned,  the  zeroth  channel  consists  of  a  direct  path 
plus  the  combination  of  all  multipath  interference  with  small  relative  de¬ 
lays  in  comparison.  The  nth  component  channel  (n  >  0)  then  represents 
multipath  interference  combinations  with  relative  delay  of  approximately 
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rn.  The  variables  .4n.  8n.  and  rn  are  modeled  as  mutually  independent 
random  vectors  with  mutually  independent  components  as  well.  That  is. 

A  =  (.4i..42. . 4i).  6  =  (6X.62 . 0L).  and  r  =  \tx.t2 . rL).  It  is  also 

assumed  that  these  quantities  are  independent  of  .40  and  8q. 

In  general,  the  range  of  the  path  delays  is  0  <  rn  <  T.  To  allow  resolvable 
demodulation  of  the  transmitted  signal,  this  rang,  is  normally  limited  such 
that  T:  <  rn  <  T  —  T-.  That  is.  the  delays  are  at  least  one  chipping  period 
awav  from  zero,  to  provide  necessary  receiver  rejection  of  the  delayed  paths, 
and  at  least  one  chipping  period  away  from  T.  so  that  the  receiver  will 
not  demodulate  the  delayed  paths  as  direct  paths  with  one-period  delay. 
However,  for  analytical  purposes,  this  paper  considers  situations  in  which 
0  <  rn  <  T  -  Tc;  the  upper  limit  is  preserved  to  model  practical  applications. 

Receiver 

The  output  of  the  channel  y(t)  is  combined  with  additive  white  Gaussian 
noise  n(t)  ~  .V  (o.  producing  the  received  signal  r(t).  The  receiver,  in 
turn,  correlates  this  signal  with  \J'2/T x(t )  cos(2~ fct)  and  integrates  from  0 
to  T.  Thus,  the  receiver  output  becomes 

Z  =  \ rz:  f  r(t)x(t)  cos('2~  fct)dt.  (6) 

V  r  j  0 

Making  the  necessary  substitutions.  (6)  becomes 


o 


Z  =  jfo  +  \/2P  b[t)x2{t)[3cos(2vfct)  +  ‘)oA0cos(2Trfct  +  $0)]cos(2iTf:1)d1 
+UT  ~nAnb(t  -  Tn)x{t  -  rn)x(t)  cos[2n  fc(t  -  Tn)  +  $n]cos(2-r  fct)dtj  | 
where  r/0  ~  A' (0,  Because  fc  is  much  larger  than  1/T  in  practical 


systems,  the  double-frequency  terms  above  may  be  ignored.  Thus,  the  inte¬ 
grated  output  becomes 


where  On  —  9n  —  '2irfcTn.  Both  6n  and  0O  are  assumed  to  be  uniformly 
distributed  in  [0.  2tt}.  Similarly,  because  of  the  size  of  fc  in  comparison  to  rn, 
the  modulo-2~  form  for  on  ~  CT[0,  2tt] .  Defining  the  multipath  interference 
component  for  the  zeroth  path  F0  =  A0bo  cos(0o)  and  splitting  the  range  of 
the  integrand  into  the  ranges  [0,Tn]  and  [ rn,T ],  (7)  becomes 

{3b0  +  70To)T  +  ~i*An  cos(°n)  (I  b(t-  rn)x{t  -  rn)x(t)dt 

n=l 

+  f  b(t  —  Tn)x(t  -  Tn)x{t)dt 

From  [3].  the  continuous  partial  autocorrelation  functions  of  x(t)  are  defined 
by 
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where  rfQ  ~  .V(0>  ^)- 


3.  Average  Probability  of  Error 


The  method  applied  in  evaluating  the  average  probability  of  error  for  the 
DS/SS  system  is  the  charact*3  istic-function  method,  as  employed  in  [1].  Here, 
the  approximation  is  based  on  the  integration  of  the  characteristic-functions 
of  the  noise  and  multipath  interference  components  of  the  correlated  receiver. 
The  analysis  which  follows  considers  deterministic  gain  coefficients  gn. 

The  gain  coefficients  are  normalized  such  that  3  =  1  and  An  does  actually 
equal  1  for  0  <  n  <  L.  From  (10),  define  the  following: 

1  =  I'Jy/Ei  (11) 


I  =  E-t  ,f».  (12) 

71=0 

This,  in  turn,  implies  that 

Z  =  y/Eb{v  -f  -f  /).  (13) 

Since  the  decision  of  the  receiver  is  based  on  comparing  the  output  Z  with  a 
zero  threshold,  the  average  probability  of  error  conditioned  on  b0  is  written 
as 

Pe  =  l-P(Z<0\bo  =  +l)  +  ^P(Z>0\bo  =  -l) 

=  l-[P(v  +  I<-l)  +  P(v  +  I>\)} 

=  \  <!?  +  /<  !)•  (14) 
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As  demonstrated  in  Appendix  A,  this  expression  becomes 


Pe  =  Q(a  1)  +  -  /  sinc(n/7r)or)(u)[l  -  0/(u)]du  (15) 

7T  JO 

where  Q(x )  =  exp(— t2/2)dt,  a2  =  (2^/A'o)"1  is  the  variance  of  7?. 

and  Or,  and  0/  are  the  respective  characteristic  functions  of  7  and  I.  Since  77 
is  a  Gaussian  random  variable,  its  characteristic  function  is  given  by 

o„(u)  =  exp  (jurj  - 

=  exp(-|uV)  (16) 

recalling  that  77  is  zero  mean.  The  characteristic  function  of  the  interference 
sum  /.  can  be  expressed  as  the  product 


©/  =  0o(u)€>(u )  (17) 

where  ©o(u)  represents  the  characteristic  function  of  the  direct-path  faded 
component  70/0  and  o(u)  denotes  the  characteristic  function-for  the  inter¬ 
ference  components  £„=1  7 nFn.  For  the  direct-path, 


0o(u)  =  E{t 

)xp(juy0Fo)] 

1 

[2t 

27T 

/  exp(j  U7060  cos  0)d$ 

Jo 

1 

t2r 

27T 

/  cos(u7o&o  cos  6)d$ 
Jo 

2 

rr/2 

=  rJ 

f  cos(7o u  cos  6)d0 

0 

where  the  sine  component  of  the  complex  exponential  integrates  to  zero. 
From.  (10),  there  are  two  cases  to  consider  in  evaluating  o(u):  6_  1  =  b0  and 
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6_!  ^  bo ,  both  of  which  are  equally  likely  to  occur.  Also,  as  opposed  to 
the  path  delay  restriction  of  [1].  r  is  assumed  to  be  uniformly  distributed  in 
[MiTc.  M2TC],  where  Mi  and  M2  are  integers  and  Mi  <  M2.  Thus. 


fa)  =  nE 


exp  ( ju  ^2  lnFn(b-i  ~  b0) 


n=l 


+  \E 


exp  ju  ^2  lnFn(b-i  #  bo) 


n=l 


=  [exP  {juinFn(b- 1  =  bo)))  +  ^  n  F  [exp  (ju“,nFn(b-i  #  bo))] 

“  n=l 


2  AA 

c  n=l 

1  L  r2-  rM2Tc  .  \ 

=  T  n  M  /  /  exP  Ju"/n  (-Rx(t)  +  ^x(r)) 

2  Jo  J M\TC  v  ' 


COS  o 


1  /*2r  [Mile  ,  .  v 

+  T  II  M  /  /  eXP  (-i?r(T)  + 

2n=i  •'°  L 

where  the  constant  M  is  defined  by 


cos  o 


drdo 

drdo 


(19) 


M  = 


1 


2t r(Ma  -  A/i)Tc 

Because  /  has  symmetric  distribution,  the  characteristic  functions  Oo  and  o 
must  be  real  and  even:  i.e..  6(u)  =  o(—u).  The  autocorrelation  functions  Rx 
and  Rx  are  defined  in  [3]  in  terms  of  the  discrete  aperiodic  autocorrelation 
function  Cx  of  the  signature  sequence  (i,).  and  the  chip  waveform  partial 
autocorrelation  functions  Rv  and  R v.  That  is. 

Rx{r)  =  Cx(N-1)R^(t-ITc)  +  Cx(X-1-1)Rv(t-ITc) 

and 


R*{t)  =  Cx{1)R*(t-ITc)  +  Cx(1  +  1)Rv(t-ITc) 
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where  /  denotes  the  integer  part  of  t/Tc.  Rk.(s)  =  /5Tc  v(t)v(t  —  $)dt.  and 
Rv(s)  =  jJc  u>(t)v(t  +  TC  —  s)dt  for  0  <  s  <  Tc.  Thus.  (19)  becomes 

o(u)  =  l-f[M  r  [M2Tc  cos {U7n[(c*(.v  -  l  -  1)  +  Cx(l  +  l))Rv(T  -  ITc) 

2  7o  JMx  tc 

a  f* OS  (5  I 

+  (Cx(Ar  - 1)  +  Cx(1))Rv{t  -  lTc))~j^\  drdo 

1  k  /*2r  rMiTc 

+  zI[M  /  cos{u7n[(-CI(A'-/-l)  +  CI(/  +  l))i?v(r-/Tc) 

2  n=l  -/°  ^MlTc 

+  (-c,(.v  -  /)  +  cx(i))RUr  -  j  drdo. 

From  [4],  the  periodic  (even)  and  odd  autocorrelation  functions  for  the  chip 
waveform  are  respectively  defined  by 

6X{1)  =  Cx(l)  +  CX{1  -  Ar),  0  <  /  <  A'  —  1 

and 

ex{l)  =  Cx(l)  -  CX{1  -  A"),  0<l<  N  -l. 

Making  these  substitutions, 


6(u)  =  l-Y[M  r  [M2Tc  cos  (u7n  \ex(l  +  1  )Rv(t  -  ITC)  +  9x(l)R^(r  -  ITC) 
2  Jo  JM\TC  l  1 

1  -hr  f  [MiTc  f  r  -  - 

+  oIlW  /  cos  U7„  \ex(l  +  l)^(r  -  1TC)  +  0x(l)R4r  -  lTe) 
Z  n=1  JO  JMiTc  (  1 

for  0  <  l  <  A’  —  2.  Considering  the  variable  substitution  s  =  t  —  ITC, 


cos  <t>  \  .  , 
~jT  f  dTd<p 

cos  o 1  .  . 
—j—  >  drdo 
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o(«)  = 
+ 


+ 


1  L  r2-  Mi-1  rTc  (  r  '1  COSO) 

in-W  £/  cos{u^n  ox(/+i)flv(S)  +  MW)j-r 

2  n=l  7°  l=M\  ”"°  t  j 

—  TT  Af  /  £  /  cos  <  U7„  !#*(/  4-  i)/M-s)  +  9x(l)RAs)\ 

2  n=l  J°  it m'°  l 


fT‘  1 

f 

/  cos  ( 

«7n 

Jo  1 

l 

jM 

|u7n 

$*(/  +  l)/?v(s)  +  @x{l)Rv{s 


T 

cos  o  | 

"T"j 

COS  o) 

Tj 

cos  ol 


1  L  Mi  —  1  r-r/2  rTc  (  »  *  .  .1  COSOl 

-  n 4M  £  /  /  cos  u7n  M + *.(/)*(-)] — } 

2  n=l  !=Mi  ^0  l  ' 


dsdo 

}  dsd© 

|  dsdo 
dsdo 


si" 


M2  ”1 


M2  -  Ml  fa 


£  f(ynu:l.Ox ) 


+  n 

n=l 


Afc-l 


M2  -  Mi  tfa 


fiw.LOx) 


where 

f(v.lh) 


Sf['2  cos  {u  [h(l  +  1  )RM  +  h(l)Rv(s)}  ^  J  dsdo.  (21) 


For  the  chip  waveform  described  in  (1),  the  partial  autocorrelation  functions 
become  Rtj{s)  =  s  and  i?v(s)  =  Tc  —  s.  0  <  s  <  Tc  [2].  So.  for  the  BPSK 
system,  the  function  in  (21)  becomes 


9  rv/2  rTc  (  ,  ,,  cosol  ,  , 

f(v:lh)  =  Zf  JQ  JQ  cos  jv^  +  1)s  +  h(l)(Tc  ~  ■»)]  — jT" |  dsdo 

2  rT/2  sin  [vh(l  +  l)T’cr~1  cos  o]  —  sin  [i>fe(/)Tcr~1  cos  o]  ^ 
“  Tfjo  v  [A(/  +1)  —  ^(03  T~l  cos  0 

Applying  the  trigonometric  identity 

sin(2A)  -  sin(2f?)  =  2  sin(.4  -  B)  cos(A  +  B), 
this  expression  reduces  to 


(20) 
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2  r/22sin{f[A((+l)-M')]^}cos{l(M/+l)  +  M')]^}J 

=  'Jo  - vw+D-mw - 

=  ;  sinc|u[)i(/+ 1)  -  ^(/)!^^|cos|f[A(/+ 1)  +  fc(/)]^?j  rfo  (22) 

which  is  the  form  obtained  in  (14)  of  [2], 

For  comparative  analysis,  the  Gaussian  approximation  of  the  probability 
of  error  is  derived  by  applying  the  average  signal-to-noise  ratio  at  the  receiver 
output  when  there  is  no  power  in  the  faded  component  of  the  zeroth  path 
(i.e..  70  =  0).  As  in  [3],  SNRq  is  defined  by 


SNRo  =  E[Z\b0  =  +l]{Var(Z|60  =  +1]}-1^2 


=  ^h[E  1  +  Xj ^Fn  |  |var  ??'  +  yfFh  ^1  +  £ fnFn j  | 


1  -1/2 


= 


Var(r?o)  +  £&Var 


r  l 


Ei,f, 


Ln=l 


! 


-1/2 


(23) 


Since  E  [eL.  7 „F„]  =  0,  Var  [Ej.,  7»F„]  =  ELi  llE(F^),  where 

E(Fl)  =  £{T-2[6_,//I(r„)  +  fl,(r„)]!cosI(®„)} 

=  ir-^dJ-.^r.)  +  fl,(T„)l2} 

=  lr-2£(fl2(r„)  +  fl2(r«)] 

-  \T-2mhaa, dZy*(T)+k'(T)]iT- 

Making  this  substitution,  (23)  becomes 
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SNR°  -  {  2  +  Eh  2T2TC(M2  -  Ml)  ln  U  ^ Rx ^  +  ^(r^r} 

_  jje + 1 _ i —  t  1 1  r2Tc\Ri(r) + ^(r)]*} 

"  \2E,  2  T2TC{M2  -  Mi)  £  ln  JmiTc  1  x{  A  J 

=  j(^)'1  +  ^[r2J,(M2  -  Mi)]-1  |>n  £>> +  **(r)1<fr}  ‘ (24) 

Substituting  the  expressions  for  i?z  and  i?x  ,  the  sum  within  the  integral 


becomes 


RI(t)  +  RI(t)  =  [Cr(.v  -  l)R*(r  -  ITC)  +  C,(.V  -  /  -  l)fl„(r  -  frt)]2 

+[C,(Z)£v(t  -  ;r£)  +  C,(i  +  l)H»(r  -  /Tt)]2 

=  kl(r  -  zrj[cj(.v  -  0  +  cj(0) 

-t  2  it  I  t  -  irc)fi*(T  -  !Tc)[Ct(.V  -  /)C','.V  -  i  -  1)  +  C,(i)C,(i  +  1)] 
+i!2  (r  -  /re)[C,2(A-  -  /  -  1)  +  Cl(l  +  1)]. 


Again,  letting  s  =  t  -ITC  and  realizing  the  expressions  for  Rv{s)  and  Rv{s), 

j M2T\R2x{t )  +  R2x(r)\dr 

J  Tc 

Mo-1 

=  Y.  {(Tc-  *)2[C,2(Ar  -  i)  +  cr2(0)  +  2  (Tc  -  S)(S)[C,(.V  -  0C.(A-  -1-1) 

l=Mi 

+C,(0C,(Z  4- 1)]  +  J2[C2(,V  -  /  -  1)  +  Cl(l  +  1)]} 

M2- 1  f  (T  -  s Y3 

=  E  -[Cx(Ar  -  0  +  ^(0]— 

(=Mi  i  ° 

/J52  53\ 

+2 [CX(N  -  l)Cx(N  -/-!)  +  CX(/)CX(Z  +  1)]  -  jj 
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+  [Cr(-V  —  /  —  !)  +  Cl([  +  1)] 


£!p 

3  i  o 


A/2-1  /  /T3\ 

=  t  2  7  [C,(.V  -  ')C-,(-V  -  /  -  1)  +  C,(1)C,(I  + 1)] 

l=Mi  l  V  0  / 

+  y  [c|(  x  -I-1)  +  Cl(i  + 1)  +  cl(.x  - 1)  +  cj(0]| . 

Thus,  the  signal-to-noise  ratio  reduces  to 


SNRo  = 


(  /9 F,\~1  1  t.  -Wj-i  j3 

=  1  V)  +  -  ^rc)]-  Y,  il  Z  tIc'(-v  -  ')c,(.v  -  /  - 1) 

t  V  ‘V°  7  2  n=l  (ssA/j 

+  CX(/)CX(/  +  1)  +  C2(A'  -  /  -  1)  +  c2x(l  +  1)  +  c2x(N  -1)  +  c2M}’1/2 

f  /2 Ek\-1  L~y2  Mi~l 

=  U^J  +  -f  E(c.(Ar-OC,(Jv-/-i)+.cI(i)cI(/+i) 

+  C*(A'  - '  -  U + <?(' + » + «?(*  -  0 + 


Ao  y 


E  [C*(.V  -  l)Cx(N-  l-  1)  +  CX(/)CX(/  +  1) 


1VAW  6N2(M2  -  Mi)  lfj^i  v 
+  c2x(N  -l-  1)  +  c2(/  + 1)  +  C2(.V  -  /)  +  Cr2(/)]}"1/2 . 

When  consideration  is  given  to  the  case  70  #  0.  the  Gaussian  approximation 
to  the  probability  of  error  must  be  conditioned  on  the  phase  shift  60,  as  seen 
in  the  expression  for  F0 ■  Thus,  as  in  [l],  the  approximation  becomes 


Pg  =  /  Q[SNRq(1  +  7Ocos0)]d0. 

Ztt  Jo 


4.  Coded  Performance  Analysis 


The  performance  of  the  spread-spectrum  system  is  enhanced  with  the 
addition  of  channel  coding,  as  one  may  anticipate.  Due  to  the  hardware 
constraints  for  physical  application  of  such  systems,  it  becomes  necessary  to 
quantize  the  continuous  receiver  output  Z  in  order  to  examine  the  effects  of 
coding.  We  start  by  defining  W  as  a  quantization  function  of  Z\  that  is, 

W  =  ivj  for  a}<Z  <  a:+1.  j  =  l,2,...,J 

where  oq  =  — oc  and  aj  =  +oc.  Here,  a.j  represents  the  lower  quantization 
threshold  for  quantization  level  j  and  we  assume  perfect  interleaving.  For 
coding  with  soft  decisions  there  are  many  levels  j ,  while  with  hard  decisions. 
j  is  only  as  high  as  two.  It  follows  from  (13)  then  that  the  probability  mass 
function  (pmf)  for  W  is  given  by 

P{W  =  Wj}  —  P{CLj  <  Z  <  Oj+i) 

=  P[aj  <  VE(tj  +  b0  +  I)  <  <Zj+1].  (27) 

Here,  Eb  is  replaced  by  E.  the  energy  per  code  symbol:  in  general,  for  rate-r 
coding,  E  =  rEb-  Letting  o  represent  the  characteristic  function  for  the  sum 
r)  +  7.  and  noting  that  <j>  is  an  even  function, 

P(a,j  <  Z  <  a,j+i\bo)  =  —f  <p(u)  cos(gu)dudg.  (28) 

J-jlm-bo  7T  Jo 

As  shown  in  Appendix  A,  one  may  replace  <p(u)  =  ov(u)oj(u)  with  Or,(u)  — 
0,,(u){l  —  d>i(u)]  to  avoid  subtracting  numbers  large  relative  to  their  difference. 
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Thus,  (28)  becomes 


P(ai  <  Z  <  aJ+1\b0) 

where  Or,  and  O/  are  given  in  (16)-(20)  and  the  function  d(u)  is  defined  by 

d(u)  =  sin  ~  t>oj  u  ~  s^n  ~  j  u  ' 

For  channel  coding,  an  upper  bound  to  the  probability  of  error  may  be  de¬ 
termined  using  D-parameter  analysis.  Denoting  the  conditional  probability- 
mass  function  byr  p(ity  |60),  the  distance  measure  is  approximated  by 

D  =  +  (30) 

j=i 

for  J  quantization  levels.  The  probability  error  bound  is  then  expressed 
as  a  polynomial  function  of  D.  For  example,  for  a  channel  with  rate- 1/2 
convolutional  code, 

Pt  <  36D10  +  21  ID12  +  1404D14  +  11633Z?16  +  77433D18  4-  502690D20  (31) 

while  for  convolutional  coding  with  rate  1/3, 

Pe  <  UD18  +  32 D20  -f  195D22  +  564D24  +  1473D26,  (32) 

where  the  weight  enumerator  polynomials  are  constructed  from  Tables  I  and 
II  of  [6]. 


=  —  f  u  1d(u)o?,(u)du 

7T  JO 

-  -/  u_1d(u)o^(u)[l  -  Qi(u)]du  (29) 

77  Jo 
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5.  Numerical  Results 


For  all  of  the  following  examples,  the  signature  sequences  were  gener¬ 
ated,  from  the  method  for  auto-optimal  (AO)  /  least-sidelobe-energy  (LSE) 
maximal-length  (m)  sequences,  as  described  in  [4].  Specifically,  the  parame¬ 
ters  listed  in  Figure  A.l  [4]  (or  [5]  when  necessary)  w'ere  applied  to  a  linear 
feedback  shift  register  to  produce  the  desired  binary  chipping  code  (i.e.,  -f  1 
and  -1  values).  These  sequences  were  chosen  for  their  structure  simplicity 
and  good  correlation  properties  [4].  As  in  [1].  it  is  assumed  that  the  relative 
power  in  the  nth  path  72  =  72  for  1  <  n  <  L.  but  the  faded  component 
power  in  the  zeroth  path.  7o.  does  not  necessarily  equal  72.  In  fact,  for  many 
instances  7^  will  assume  a  zero  value  (i.e..  no  faded  component  in  the  ze¬ 
roth  path).  All  of  the  numerical  results  were  evaluated  using  the  computer 
programs  listed  in  Appendix  B. 

The  first  example  considers  the  effect  on  the  average  probability  of  er¬ 
ror,  determined  from  the  characteristic-function  method  (15).  when  the  path 
delay  is  less  than  a  chipping  period.  Figure  1  shows  a  comparison  in  er¬ 
ror  probability  between  two  PSK  DS/SS  systems,  each  multipath  channel 
with  L  =  5  reflected  paths,  individually  having  relative  power  7 2  =  0.4.  and 
7q  =0.1.  The  signature  sequence  used  is  the  X  =  31  sequence  taken  from  the 
first  line  of  Figure  A.l  [4].  In  the  first  system,  the  path  delay  r  is  uniformly 
distributed  in  [Tc,  T  —  Tc],  while  in  the  second,  the  delay  ranges  uniformly  in 

[o.r-Tc], 
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According  to  the  figure,  the  system  performs  well  when  the  path  delay 
is  at  least  one  chipping  period:  the  probability  of  error  quickly  decreases 
for  increasing  E^j. Yo  values.  However,  when  r  ~  L  [0.  T  —  77] .  the  error 
probability  decreases  at  a  much  slower  rate  and  eventually  levels  off  at  about 
1CT4.  This  demonstrates  the  fact  that  in  the  second  system  the  multipath 
interference  component  arrives  too  close  to  the  zeroth  path  signal.  As  a 
result,  the  receiver  is  unable  to  reject  the  delayed  path  effects. 

The  second  example  is  a  comparison  between  narrowband  and  broad¬ 
band  spread-spectrum  systems.  In  the  first  case  (narrowband),  a  signature 
sequence  of  length  .V  =  127  is  created  using  shift  register  polynomial  param¬ 
eter  =  301  and  initial  register  loading  a0  =  0010010  [ok  The  multipath 
channel  has  5  delayed  paths,  each  with  r  ~~  7’[0. 1077] .  For  the  broadband 
case,  a  sequence  of  length  .V  =  1023  is  generated  using  H  =  2201  and 
q0  =  0111011010  [5].  Here,  each  of  the  5  paths  is  uniformly  distributed  in 
[0.8077],  This  corresponds  to  the  same  delay  spread  if  the  data  rate  is  as  in 
the  .V  =  127  case.  In  both  cases,  the  faded  component  of  the  zeroth  path  is 
removed  (i.e..  *,q  =  0).  and  each  path  has  relative  power  *-2  =  0.4.  Figure  2 
shows  the  probability  of  error  for  each  application. 

As  one  might  expect,  the  broadband  application  outperforms  the  narrow- 
band  one  for  higher  values  of  £<,/ A’o-  Clearly,  as  A'  increases,  the  chipping 
period  T 7  must  decrease,  assuming  constant  data  pulse  length  T :  thus,  provid¬ 
ing  greater  multipath  rejection  capability  for  the  receiver  and.  hence,  lower 
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error  probability.  However,  such  performance  comes  at  the  expense  of  in¬ 
creased  bandwidth  and  system  complexity  (more  shift  registers  are  required 
for  sequence  generation). 

In  the  next  example,  the  average  probability  of  error  is  analyzed  for  the 
two  systems  of  Figure  2  as  a  function  of  the  number  of  reflected,  or  delayed, 
paths  £.  For  a  fixed  Eb/N0  value  of  10  dB,  Figure  3  demonstrates  that  this 
measure  of  performance  increases  as  more  paths  are  added  to  the  channel. 

For  the  fourth  example,  the  narrowband  and  broadband  spread-spectrum 
systems  are  examined  for  a  particular  multipath  environment.  There  are  a 
total  of  L  =  12  delayed  paths,  divided  into  four  groups  of  three  paths,  each 
group  with  a  particular  relative  power  for  the  individual  paths,  and  =  0. 
For  the  narrowband  (.V  =  127  chips/data  bit),  there  are: 

3  paths:  r  ~  U[0ATc\,  72  =  1 

3  paths:  r  ~  £7[24TC, 287’c],  72  =  0.1 

3  paths:  r  ~  f/[48Tc,  52rcj.  72  =  0.1 

3  paths:  r  ~  £/[96Tc,  100rc],  72  =  0.03 

For  the  broadband  (N  =  1023  chips/data  bit),  the  length  of  each  data  pulse 
T  is  maintained  from  the  narrowband  case.  Thus,  the  12  paths  for  this 
system  are  divided  as  follows: 


3  paths: 

t  ~  £/[0, 33TC], 

72  =  1 

3  paths: 

r~  t/[194rc,226Tc], 

-> 

II 

o 

) — * 

3  paths: 

r~  £'[387rc,419rc], 

72  =  0.1 

3  paths: 

r  ~  U[774Tc,8Q6Tc\. 

7 2  =  0.03 
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No.  of  delayed  paths  L 


Figure  3:  Probability  of  error  for  increasing  number  of  delayed  paths 
(Eb/X0  =  10  dB,  72  =  0.4,  7q  =  0). 


Error  Probability,  Pe 


Figure  4:  Probability  of  error  for  a  particular  multipath  environment  ( L  =  12 
total  paths,  7o2  =  0). 
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As  shown  by  Figure  4,  while  the  probabilities  of  error  for  both  cases  are 
fairly  close  at  low  Eb/ N0  values  (0-3  dB),  their  relative  difference  does  in¬ 
crease  for  higher  Eb/N0.  Notice,  though,  that  these  systems  perform  rather 
poorly  in  comparison  to  previous  situations;  Pe  lies  above  10“4  even  at 
Eb  /No  =  16dB.  This  is  due  to  several  factors.  First,  the  path  delay  re¬ 
gion  in  which  a  portion  of  multipath  rejection  is  unattainable  (i.e.,  r  <  Tc) 
also  has  the  greatest  amount  of  relative  power  for  each  path  (72  =  1).  As 
a  result,  this  interference  component  receives  the  most  emphasis  of  the  four. 
Second,  the  high  number  of  paths  in  total  (12)  should  generally  imply  higher 
error  probability,  as  proven  from  Figure  3. 

Figure  5  compares  three  different  systems  for  the  length  N  —  31  signature 
sequence  in  which  path  delay  is  uniformly  distributed  in  [0,  Tc],  =  0,  and 
the  sum  of  the  multipath  relative  powers  is  constant  for  each  system  (i.e., 
J2n=\  7n  =  Tl2  =  0-2)-  The  total  number  of  paths  L  increases  from  2  to  100, 
as  shown.  The  L  =  100  curve  coincidesAthe  Gaussian  approximation  to  the 
error  probability,  expressed  in  (26),  which  depends  only  on  7«- 

The  impact  of  the  multipath  interference  is  demonstrated  by  the  signif¬ 
icant  difference  between  the  L  =  2  data  and  that  of  L  =  10  and  100  for 
Eb / A o  >11  dB.  More  importantly,  the  figure  shows  that  the  performance  of 
the  characteristic-function  method  approaches  that  of  the  Gaussian  model 
for  high  values  of  L.  In  the  case  of  L  —  10,  the  characteristic-function  model 
appears  Gaussian  up  to  about  Eb  /No  =  13  dB.  while  for  L  =  100.  the  ab- 
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Figure  5.  Probability  of  error  for  characteristic-function  method  approaching 
Gaussian  approximation  in  performance  (N  =  31,  r  ~  U[0.Tc}.^  =  0). 
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solute  difference  in  error  probability  is  less  than  2  x  10-8.  From  the  Central 
Limit  Theorem,  one  would  expect  that  for  sufficiently  high  L  the  distribu¬ 
tion  of  the  multipath  interference  is  indeed  approximated  by  the  Gaussian 
distribution. 

The  last  two  examples  examine  the  performance  of  DS/SS  systems  with 
and  without  channel  coding.  In  the  first,  rate-1/2  and  1/3  convolutional 
coding  is  applied  to  a  system  with  N  =  31  chips/data  bit  signature  sequence 
and  path  delay  r  ~  U[TC ,  T  —  Tc]\  thus,  the  demodulating  sequence  has  length 
N  =  62  or  93  chips/  coded  symbol,  depending  upon  the  chosen  rate.  Figure  6 
compares  the  probability  of  error  bounds  (31,32)  for  the  coded  systems  with 
the  average  error  probabilities  of  uncoded  systems  with  effectively  similar 
chipping  sequence  lengths.  That  is,  the  rate-1/2  system  is  plotted  against  an 
uncoded  one  having  sequence  length  N  =  63  chips/data  bit  (H  =  103  and 
a0  =  000010  [4])  and  r  ~  t/[3Tc,  61TC],  while  the  counterpart  for  the  rate- 
1/3  system  has  N  =  127  chips/data  bit  and  r  ~  U[5TC,  1232^] ,  assuming 
the  same  data  bit  duration  T  of  the  original  system.  Both  cases  of  soft 
and  hard  decisions  are  included,  using  (31)-(32).  For  soft  decisions,  a  large 
number  (>100)  of  quantization  levels  was  chosen.  The  figure  also  shows  the 
error  probability  for  the  system  prior  to  any  channel  coding  as  a  matter  of 
reference.  Comparisons  of  coded  and  uncoded  systems  is  somewhat  biased 
because  the  performance  of  the  coded  system  is  based  on  an  upper  bound, 
while  the  un coded  system  is  based  on  an  exact  calculation. 


27 


Figure  6:  Probability  of  error  for  DS/SS  systems,  with  and  without  channel 
coding  (L  =  5,72  =  0.4, -/q  =  0). 


As  expected,  the  coded  systems  achieve  lower  error  probabilities  for  high 
Eb/X o  than  those  of  the  comparison  models.  However,  at  Eb /X0  less  than 
about  5.5  dB.  the  rate-1/2  coded  system  with  hard  decisions  starts  to  perform 
worse.  When  Eb/X0  is  as  low  as  about  1.2  dB.  even  soft  derisions  with  the 
rate-1/3  coded  system  fail  to  perform  effectively.  This  is  due  to  the  fact  that 
at  sufficiently  low  input  signal-to-noise  ratios,  channel  coding  exploits  the 
weak  energy  of  the  data  signal  and  introduces  further  errors  to  the  system. 
As  a  result,  it  is  not  uncommon  to  even  find  error  probabilities  greater  than 
one  (e.g.,  rate-1/2  coding  (hard  decision)  at  Eb/X0  =  0  dB).  For  the  most 
part,  the  systems  with  soft  decisions  perform  about  3  dB  better  than  hard 
decision  schemes. 

Finally,  the  channel  coding  schemes  above  are  applied  to  the  fourth  ex¬ 
ample  situation  in  which  there  is  a  particular  multipath  environment  for  each 
system.  Figures  7  and  8  show  the  same  type  of  probability  of  error  data  as 
in  the  previous  example  for  the  narrowband  and  broadband  spread-spectrum 
systems,  respectively. 

For  the  narrowband  case  (.V  =  127).  an  uncoded  system  with  chipping 
sequence  length  N  =  2bb(H  =  455  and  q0  =  01101111  [4])  is  compared  with 
a  rate-1/2  coded  system  (N  =  254  chips/information  bit).  Coding  with  soft 
decisions  perform  better  after  Eb/Xo  =  3.8  dB.  wffiile  hard  decisions  pro¬ 
duce  lower  error  probability  only  after  Eb/No  =  7.3  dB.  The  rate-1/3  coded 
systems  (X  =  381  chips/information  bit)  are  analyzed  against  an  un coded 
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Probability  of  Error 


Figure  8:  Probability  of  error  for  broadband  system  and  particular  multipath 
environment,  with  and  without  channel  coding  ( L  =  12  total  paths,  7q  =  0). 
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system  with  signature  sequence  X  =  oll^H  =  1041  and  Qo  =  010101001 
[5]).  Here,  for  Eb/Xo  >1.7  dB,  the  soft-decision  coding  outperforms  the 
uncoded  system,  while  hard  decisions  take  effect  only  after  Eb/Xo  =  4.7  dB. 
The  threshold  distinguishing  the  effectiveness  of  soft  decisions  over  hard  de¬ 
cisions  is  at  Eb/Xo  —  2.4  dB.  prior  to  which  the  rate-1/3,  hard-decision  code 
does  better  than  that  with  rate-1/2  soft  decisions. 

For  the  broadband  case  (X  =  1023),  the  rate-1/2  coded  system  (A*  = 
2046  chips/information  bit)  is  compared  with  an  uncoded  system  with  chip¬ 
ping  sequence  length  X  =  2047 (H  =  5001  and  oq  =  11000001000  [5]).  At 
Eb/Xo  =  2.8  dB,  soft  decisions  begin  to  produce  lower  error  probability,  while 
the  same  is  true  for  hard  decisions  at  Eb/Xo  >  5.5  dB.  The  rate-1/3  coded 
system  has  bandwidth  comparable  to  an  uncoded  system  with  X  =  4095 (H 
=  14501  and  a0  =  101111000001).  In  this  case,  for  Eb/Xo  >  1.2  dB.  the 
coded  system  with  soft  decisions  performs  better  than  the  uncoded  compar¬ 
ison,  while  hard  decisions  do  not  take  effect  until  after  Eb/Xo  =  3.9  dB.  For 
this  particular  system,  the  rate-1/3  code  with  hard  decisions  is  more  effective 
than  the  rate-1/2  code  with  soft  decisions  for  Eb/Xo  <  1.1  dB. 
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6.  Conclusions 


From  Figure  1.  it  is  clear  that  the  DS/SS  system  is  only  effective  at  reject¬ 
ing  multipath  when  the  path  delays  are  at  least  one  chip  period.  One  method 
of  reducing  the  average  probability  of  error  is  to  increase  the  processing  gain 
X.  However,  this  comes  at  the  expense  of  bandwidth  and  transmitter  com¬ 
plexity  (i.e.,  more  shift  registers  required  for  signature  sequence  generation). 
When  more  delayed  paths  are  added  to  the  system,  the  multipath  interference 
increases  and  system  performance  degrades  rapidly.  For  a  sufficiently  high 
number  of  paths,  the  performance  obtained  via  the  characteristic-function 
method  approaches  that  of  the  Gaussian  approximation,  as  expected  from 
the  Central  Limit  Theorem. 

The  DS/SS  system  is  significantly  enhanced  with  the  application  of  chan¬ 
nel  coding.  In  fact,  according  to  Figure  6.  the  probability  of  error  drops  to 
lO-10  at  about  Eb/Xo  =  6.3  dB  for  a  system  with  rate- 1/3 -convolutional 
coding  and  soft  decisions  applied.  For  the  spread-spectrum  sysiems  with  a 
particular  multipath  environment,  there  is  a  significant  difference  between 
the  narrowband  and  broadband  cases  at  low  error  probability  (high  Eb/Xo). 
However,  for  relatively  higher  probability  of  error  (low  Eb/Xo),  the  broad¬ 
band  system  provides  a  1-2  dB  improvement  over  the  narrowband  counter¬ 
part. 
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Appendix  A 

Characteristic-Function  Approach  To  Probability  of  Error 

In  this  appendix,  for  purposes  of  completeness,  we  derive  the  expression 
in  [1]  for  the  average  probability  of  error,  given  in  terms  of  the  characteristic 
functions  of  the  multipath  interference  and  noise.  From  (14)  of  the  text,  the 
error  probability  is  expressed  as 

P.  =  i-  |f(-l<ntl<l).  (33) 

Letting  G  =  rj  +  I ,  and  denoting  the  probability  density  of  G  by  /c(p), 
(33)  becomes 

P,  =  \  -  \  J\  fc(9)<l9.  (34) 

Since  the  random  variables  tj  and  I  are  mutually  independent  and  have 
symmetric  distributions,  their  characteristic  functions,  as  well  as  that  of  their 
sum,  must  be  real  and  even.  Denoting  the  characteristic  functions  of  tj.  I. 
and  tj  +  /  by  <pv,  <z>j,  and  <p,  respectively,  fdg)  is  given  by 

/g(s)  =  ^  j_^e~3u96{y)du 

=  —  c>(u)cos  (gu)du 
T  Jo 

since  fdg)  must  also  be  real-valued.  Substituting  into  (34), 

1  1  r1 2  1 

Pe  =  t-x/  ~  Q(u)  cos{gu)dudg 

2  2  J- 1  7T  Jo 
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1  1  roc  r\ 

= - /  /  o(u)  cos(gu)dgdu 

2  Jo  Jo 

=  -  —  —  f  o(w)u-1  sin(i/)du. 

2  TT  Jo 

To  avoid  subtracting  very  large  numbers  relative  to  their  difference,  one  may 
replace  <t>{u )  =  4>r,Oi  with  0n  —  —  <Pi}-  Thus, 

1  1 

Pe  =  - /  tT1  sin  (u)<pn{u)du 

1  IT  JO 
1  roo 

+  -  u_1  sin(u)0,,(u)[l  -  0j(u)]du.  (35) 

7T  Jo 

In  the  absence  of  interference  (i.e.,  7  =  0),  0/  =  1  and  error  is  due  only  to 
the  presence  of  Gaussian  noise.  Denoting  this  noise  error  probability  by  Pn. 
(34)  becomes 

pe  =  pv  +  jt-1  /  tT1  sin(u)0t7(u)[l  -  <pi(u)\du.  (36) 

Jo 

Since  rj  is  Gaussian,  Pv  is  given  by  Q(cr_1),  where  a2  is  defined  as  the 
variance  of  77,  The  result  is  (15)  of  the  text. 
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Appendix  B 

Computer  Programs  For  Characteristic-Function  Method 

The  following  programs  provide  simulations  for  the  average  probability 
of  error  (15),  the  Gaussian  approximation  (26),  and  error  bounds  for  cod¬ 
ing  (31,32).  The  computer  code  is  written  in  MATLAB  version  4.0  and  all 
integrations  are  approximated  applying  Simpson's  method. 

The  first  program,  PROBERROR.m,  computes  the  average  probabil¬ 
ity  of  error  for  both  the  characteristic-function  method  and  the  Gaussian  ap¬ 
proximation.  The  user  is  able  to  select  the  desired  signature  sequence  from  a 
menu;  the  program  assumes  existing  data  files  containing  the  sequence  itself 
and  all  the  appropriate  autocorrelation  functions.  The  user  may  also  create 
one  or  more  multipath  environments,  each  with  a  desired  range  for  the  path 
delay.  The  following  function  files  are  called  within  the  main  program: 

SIMP.m  -  A  Simpson’s  approximation  to  the  <p0  integral. 

PHI.m  -  Computes  <t>;  calls  function  file  FF.m  which  provides  Simpson’s 
method  for  function  f[y\  L  h ). 

The  second  program,  PROBMASS.m,  first  determines  the  probability 
mass  function  for  the  quantized  receiver  output  (29).  This  is  then  used  to 
compute  the  distance  parameter  D  (30),  which  is  applied  in  approximating 
upper  probability  bounds  for  channel  coding  purposes;  the  program  considers 


convolutional  coding  with  rates  1/2  and  1/3.  The  pmf  data  is  calculated  from 
the  same  two  function  files  of  the  previous  program. 
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%  PROBERROR.m 
% 

%  A  program  to  compute  the  average  probability  of  error  of  a  DS/SS  system  with 
%  multipath  fading  and  BPSK  modulation,  using  both  the  characteristic-function 
%  method  and  the  Gaussian  approximation.  Path  delay  is  uniform  from  one 
%  multiple  of  Tc  to  another  for  a  given  number  of  paths  within  each  multipath 
%  environment.  User  selects  desired  chipping  sequence  from  a  menu. 

% 

%  Input:  gamma0''2  (assumed  constant);  for  each  path  delay  situation  —  number  of 
%  paths  L,  gamma~2,  uniform  range  for  path  delay;  Eb/No,  upper  limit  and 

%  number  of  points  of  integration  for  Simpson's  method. 

% 

%  Output:  for  each  Eb/No,  average  probability  of  error  for  both  characteristic- 
%  function  method  and  Gaussian  approximation. 

% 


clear;  format  long  e; 

k  =  menu ( 'Choose  sequence  length ', 'N=31 ', 'N=63 ', 'N=127 ', 'N=255 ', 'N=511 ',.. . 
'N=1023 ' , 'N=2047 ' , 'N=4095 ' ) ; 


if  k  ==  1, 

load  code_31 

elseif  k  ==  2, 

load  code_63 
elseif  k  ==  3, 

load  code_127 
elseif  k  ==  4, 

load  code_255 
elseif  k  ==  5, 

load  code_511 
elseif  k  ==  6, 

load  code_1023 
elseif  k  ==  7, 

load  code_2047 

else 


%  each  signature  sequence  file  contains  the  actual 
%  sequence  (code)  and  all  autocorrelation  functions: 
%  aperiodic:  R  periodic ( even ) :  theta 

%  odd:  thetahat 


load  code_4095 

end 

disp( 'Enter  the  requested  data:') 

gammaOsq  =  input ( 'gamma_0-squared  =  ');  %  assumed  to  be  constant  for  all 

%  multipath  environments. 

gammaO  =  sqrt (gammaOsq) ; 

sit  =  input ( 'How  many  different  delay  situations  are  there?  '); 
for  j  =  l:sit 

disp( 'situation  #');  disp(j); 

L ( j )  =  input ('no.  of  paths  L  =  '); 
gammasq(j)  =  input ( 'gamma- squared  =  '); 
disp( 'Path  delay  tau  -  U[Ml*Tc,  M2*Tc]:') 

Ml ( j )  =  input ( 'constant  multiple  Ml  =  '); 

M2(j)  =  input ( 'constant  multiple  M2  =  '); 
gamma(j)  =  sqrt(gammasq( j ) ) ; 

end; 

N  =  length ( code ) ;  %  length  of  sequence 


data  =  input ( 'How  many  Eb/No  data  values?  '); 
ssuml  =  zeros ( 1 ,data) ;  ssum2  =  zeros ( 1 , data) ; 
for  d  =  l:data 

EbN0_dB(d)  =  input ('Enter  Eb/NO  value  (dB):'); 
b  =  input ( 'upper-limit  of  integration  =  '); 
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n  =  input('no.  of  pts.  (must  be  even)  for  integration  =  '); 
h  =  b/n;  EbNO(d)  =  10" (EbN0_dB(d)/10 ) ; 

sigmasq(d)  =  l/( 2*EbN0 (d) ) ;  sigmainv(d)  =  1/sqrt ( sigmasq(d) ) ; 

Q(d)  =  Qfunct ( sigmainv(d) ) ;  %  compute  the  Q-f unction  of  sigma- 

%  inverse. 

i  =  [ 0 : n ] ;  %  Simpson's  method  starts  here,  done  using  vector  form, 

u  =  IE-100  +  i*h; 

phi_eta=  exp(-(u."2)  *  (sigmasq(d) )/2 ) ; 

phi_tilde_0  =  2/pi  *  simp(gammaO*u,  pi/2,  20);  %  SIMP.m  is  a  function 

%  which  performs 
%  Simpson's  method  for 
%  phi_tilde_0. 

phi_tilde  =  phi(gamma'*u,L,Ml,M2,N,theta,thetahat) ;  %  see  separate 

%  listing  for 
%  function  PHI.m. 

phi_I  =  phi_tilde  0  .*  phi_tilde; 

g  =  ones(l,n+l) ,/u  .*  sin(u)  .*  (l-phi_I)  .*  phi_eta; 

coeff  =  2  *  (rem( i, 2 )+l ) ;  %  coefficients  for  Simpson's  summation. 

coeff(l)  =  1;  coeff (n+1)  =  1; 
ssum(d)  =  sum(coeff  .  *  g) ; 
simpson(d)  =  h/3  *  ssum(d); 

proberr(d)  =  Q(d)  +  1/pi  *  simpson(d);  %  this  is  eq.  (15)  of  text. 

%  This  section  computes  the  Gaussian  approximation 
sum2  =  0; 
for  j  =  l:sit 

suml  =  0; 

for  lg  =  Ml ( j ) :M2 ( j ) — 1 

suml  =  suml  +  R( 2*N+l-lg) *R( 2*N-lg)  +  R(N+l+lg) *R(N+2+lg) . . . 

+  R( 2*N+l-lg) "2  +  R( 2*N-lg) "2  +  R(N+l+lg)"2  +  ... 
R(N+2+lg) "2; 

end 

sum2  =  sum2  +  1/6  *  L( j ) *gammasq( j )/(N~2* (M2 ( j )-Ml ( j ) ) )  *  suml; 

end 

SNR_0  =  1  /  sqrt(sigmasq(d)  +  sum2); 

ng  =  500;  %  choose  500  points  for  Simpson's  method  integration. 

ig  =  [0:ng];  hg  =  pi/ng;  %  may  use  pi  as  an  upper  limit  since  function 

%  for  integration  is  even 

thetag  =  ig*hg; 

g  =  Qfunct (SNR_0*(1  +  gammaO * cos (thetag) )) ; 
coeff g  =  2  *  (rem(ig,2)+l) ; 
coeffg(l)  =  1;  coeffg(ng+l)  =  1; 
ssumg(d)  =  sum(coeffg  .*  g); 
simpsong(d)  =  hg/3  *  ssumg(d); 

proberrg(d)  =  1/pi  *  simpsong(d);  %  1/pi  factor  instead  of  l/(2pi) 

%  since  integrating  only  to  pi;  this 
%  is  eq.  (26) . 

disp('ave.  prob.  of  error  for  char-funct  method  =  ');  disp(proberr (d) ) 
disp( 'Gaussian  approximation  to  prob.  of  error  =  ');  disp(proberrg(d) ) 

end 
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function  simpson  =  simp( const, b,n) 

%  SIMP  Simpson's  rule  for  phi_tilde_0  characteristic  function. 

% 

%  Input:  gamma_0*u  (const),  upper-limit  (b)  and  number  of  points  for 

%  integration  (n). 

% 

%  Output:  Simpson's  approximation  to  the  integral. 

h  =  b/n;  ssum  =  0; 
i  =  [  0 :  n  ] ; 
theta  =  i*h; 

g  =  cos (const'  *  cos (theta));  %  function  for  inegration. 

coeff  =  2  *  (rem(i,2)+l) ;  %  coefficients  for  approximation. 

coeff(l)  =  1;  coeff (n+1)  =  1; 
ssum  =  coeff  *  g'; 

simpson  =  h/3  *  ssum;  %  Simpson's  approximation 
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function  out  =  phi(var, L,M1,M2,N, theta, thetahat) 


%  PHI  Characteristic  function  phi_tilde(u) .  Function  f  for  theta  and  thetahat 
%  is  called  by  FF.m. 

% 

%  Input:  variable  (var),  no.  of  reflected  paths  (L) ,  constant  multiples  of 

%  Tc  (Ml ,M2 ) ,  length  of  signature  sequence  (N) ,  and  autocorrelation 

%  vectors  theta  and  thetahat. 

% 

%  Output:  the  value  of  the  characteristic  function. 

s  =  length (L);  %  number  of  multipath  environments. 

prodl  =  1;  prod2  =  1; 
for  j  =  l:s 

ssuml  =  0;  ssum2  =  0; 
for  1  =  Ml ( j ) :M2 ( j )-l 

f  =  ff (var( j , : ) , l,theta,N);  %  compute  f  for  theta  and 

%  thetahat 

f_hat  =  ff (var( j; : ) , 1, thetahat, N); 


\ 


end 

out 


ssuml  =  ssuml  +  f; 
ssum2  =  ssum2  +  f__hat; 

end 

outl  =  ( 1/(M2 ( i )-Ml ( j ) )  *  ssuml)  . ~  L(j); 
out2  =  ( 1/ ( M2 ( j ) -Ml ( j ) )  *  ssum2 )  . *  L(j); 
prodl  =  prodl  .  *  outl; 
prod2  =  prod2  . *  out2 ; 

0.5  *  (prodl  +  prod2);  %  this  corresponds  to  phi_tilde(u) . 


function  f  =  ff(v,k,t,N) 

%  FF  Function  f  for  vectors  theta  and  thetahat.  Considers  zeroeth  index  for 
%  autocorrelation  vectors.  Input  is  real  variable  (v),  index  (k), 

%  autocorrelation  vector  (t),  and  length  of  code  sequence  (N).  Output  is 

%  Simpson's  formula  approximation. 

b  =  pi/2;  n  =  20;  %  only  20  points  are  necessary  for  good  approximation, 

h  =  b/n; 

i  =  [0:n];  1  =  k+1; 
theta  =  i*h; 

gl  =  sinc(v'*(t(l+l)-t(l) )*cos(theta)/(2*pi*N) );  %  SINC.m  simply  performs  the 

%  sine  function  for  given 
%  variable. 

g2  =  cos(v'*(t(l+l)+t(l) )*cos(theta)/(2*N) ); 

g  =  gl  .*  g2; 

coeff  =  2  *  (rem(i, 2)+l) ; 

coeff(l)  =  1;  coeff (n+1)  =  1; 

ssum  =  coeff  *  g'; 

simpson  =  h/3  *  ssum;  f  =  2/pi  *  simpson;  %  this  is  the  function  f. 
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k  PROBMASS.m 
k 

k  A  program  to  compute  the  probability  mass  function  of  Z,  P(a<Z<b|bO),  of  a 
k  multipath  DS/SS  system  with  BPSK  modulation.  Path  delay  is  uniform  from  one 
k  multiple  of  Tc  to  another,  for  a  given  number  of  multipath  environments.  User 
k  selects  desired  chipping  sequence  from  a  menu.  Here,  gammaO-squared(gammaOsq) 
k  is  assumed  to  be  zero,  but  may  be  reset  to  constant  for  all  delay  situations, 
k  The  pmf  is  used,  in  turn,  to  determine  the  Hamming  D-parameter  and  apply 
k  toward  computing  the  probability  error  bounds  for  specific  channel  coding, 
k 

k  Input:  for  each  delay  invironment  —  number  of  paths  (L),  gamma-squared 
k  (gammasq),  Tc  multiples  (M1,M2); 

k  total  number  of  ranges  for  receiver  output  Z,  epsilon  increment,  Eb/No, 

k  upper  limit  and  no.  of  points  for  integration  using  Simpson's  method, 

k 

k  Output:  probability  error  bounds  for  rate-1/2  and  1/3  convolutional  codes 
zlear;  format  long  e; 

ft  =  menu ( 'Choose  sequence  length ', 'N=31 'N=127 ', 'N=1023 ' ) 
if  k  ==  1, 

load  code_31 
elseif  k  ==  2, 

load  code_127 

else 

load  code_1023 

end 

disp( 'Enter  the  requested  data:') 

sit  =  input('How  many  different  delay  situations  are  there?  '); 

for  j  =  l:sit  %  same  data  requested  as  in  PROBERROR.m 

disp( ' situation  #');  disp(j); 

L ( j )  =  input ('no.  of  paths  L  =  '); 
gammasq(j)  =  input (' gamma-squared  =  '); 

disp( 'Path  delay  tau  -  U[Ml*Tc,  M2*Tc]:') 

Ml ( j )  =  input ( 'constant  multiple  Ml  =  '); 

M2(j)  =  input ( 'constant  multiple  M2  =  '); 

gamma0sq=0;  %  assume  gammaO-squared  is  zero  throughout. 

gamma(j)  =  sqrt ( gammasq ( j ) ) ;  gammaO  =  sqrt ( gammaOsq) ; 

end; 

N  =  length ( code ) ; 

num  =  input ( 'How  many  Eb/No  values?'); 
for  g  =  l:num 

EbN0_dB(g)  =  input ('Enter  Eb/NO  (dB):'); 

eps  =  input ( 'Enter  increment  epsilon:');  %  increment  used  in  mapping 

%  out  the  discrete  Z  values. 

data  =  input('How  many  total  ranges  for  Z  (even)?  '); 

p  =  [-data/2 : data/2 ] ;  %  split  the  desired  ranges  into  positive  and 

%  negative  regions. 

Eb(g)  =  10-(EbN0_dB(g)/10);  N0(g)  =1;  %  the  energy  per  information  bit 

%  Eb,  is  needed  specif ically . 
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E2(g)  =  Eb(g)/2;  E3(g)  =  Eb(g)/3;  %  this  is  the  energy  per  coded 

%  symbol  for  rates-1/2  and  1/3 
%  codes . 

EN02 (g)  =  (Eb(g)/N0(g) )/2;  EN03(g)  =  (Eb(g)/N0(g) )/3; 
bb  =  input ( 'upper-limit  of  integration  =  '); 
n  =  input ('no.  of  pts.  (must  be  even)  for  integration  =  '); 
h  =  bb/n; 

sigmasq2 (g)  =  l/( 2*EN02 (g) ) ;  sigmasq3(g)  =  l/(2*EN03(g) ); 
sigmainv2(g)  =  l/sqrt(sigmasq2(g) ) ;  sigmainv3(g)  =  1/sqrt ( sigmasq3 (g) ) ; 

d2  =  eps*p*sqrt ( E2 ( g) ) ;  d3  =  eps*p*sqrt(E3(g) ) ;  %  these  are  the  actual 

%  values  for  Z . 

a2  =  d2(l:deta);  a3  =  d3(l:data);  %  a  and  b  the  range  limits  for 

b2  =  d2 ( 2 :data+l ) ;  b3  =  d3 ( 2 : data+1 ) ;  %  the  Z  values.  Throughout,  a 

%  2-suffix  refers  to  a  variable 
%  for  rate-1/2  code,  while  a 
%  3-suffix  refers  to  that  for 
%  rate-1/3  code. 

i  =  [  0  :  n  ] ; 
u  =  IE-100  +  i*h; 

uinv  =  ones(l,n+l) ./u;  uinvd  =  uinv; 
phi_eta2  =  exp(-(u.~2)  *  sigmasq2 (g)/2 ) ; 
phi_eta3  =  exp(-(u./'2)  *  sigmasq3  ( g)/2  ) ; 
phi_etad2  =  phi  eta2;  phi_etad3  =  phi_eta3; 
phi_tilde_0  =  2~/pi  *  simp(gammaO*u,  pi/2,  20); 
phi_tilde  =  phi ( gamma '*u,L, Ml, M2 ,N, theta, thetahat ) ; 
phi_I  =  phi_tilde_0  .*  phi_tilde;  phi_Id  =  phi_I; 

for  m  =  l:data-l  %  generates  data-length  column  vectors  for 

%  u-inverse,  phi_eta,  and  phi_I . 

uinvd  =  [uinvd;  uinv]; 
phi_etad2  =  [phi_etad2;  phi_eta2 ] ; 
phi_etad3  =  [phi_etad3;  phi_eta3]; 
phi_Id  =  [phi_Id;  phi_I]; 

end 

%  computes  functions  for  Simpson's  method,  considering  both  b0=+l  (pos) 
%  and  b0=-l  (neg) . 

sindiff2pos  =  sin( (b2'/sqrt(E2(g) )-l)*u)  -  sin( (a2 '/sqrt (E2 (g) )-l ) *u) ; 

sindiff3pos  =  sin( (b3'/sqrt(E3(g) )-l)*u)  -  sin( (a3 '/sqrt (E3 (g) )-l ) *u) ; 

sindiff2neg  =  sin( (b2 '/sqrt (E2 (g) )+l ) *u)  -  sin( (a2 '/sqrt (E2 (g) )+l ) *u) ; 

sindiff3neg  =  sin( (b3 '/sqrt (E3(g) )+l )*u)  -  sin( (a3 '/sqrt (E3 (g) )+i ) *u) ; 

gl2pos  =  uinvd  .*  sindiff2pos  . *  phi_etad2; 

gl3pos  =  uinvd  . *  sindiff3pos  .*  phi_etad3; 

g22pos  =  uinvd  .*  sindiff2pos  . *  (l-phi_Id)  .*  phi_etad2; 

g23pos  =  uinvd  .*  sindiff3pos  . *  (l-phi_Id)  .*  phi_etad3; 

gl2neg  =  uinvd  .*  sindiff2neg  .*  phi_etad2; 

gl3neg  =  uinvd  .*  sindiff3neg  . *  phi_etad3; 

g22neg  =  uinvd  .*  sindiff2neg  .*  (i-phi_Id)  .*  phi_etad2; 

g23neg  =  uinvd  .*  sindiff3neg  . *  (l-phi_Id)  .*  phi_etad3; 

coeff  =  2  *  (rem( i, 2 )+l ) ; 

coeff(l)  =  1;  coeff (n+1)  =  1; 

ssuml2pos  =  coeff  *  gl2pos';  ssuml3pos  =  coeff  *  gl3pos'; 
ssum22pos  =  coeff  *  g22pos';  ssum23pos  =  coeff  *  g23pos'; 


ssuml2neg  =  coeff  *  gl2neg';  ssuml3neg  =  coeff  *  gl3neg'; 
ssum22neg  =  coeff  *  g22neg';  ssum23neg  =  coeff  *  g23neg'; 
simpsonl2pos  =  h/3  *  ssuml2pos;  simpsonl3pos  =  h/3  *  ssuml3pos; 

simpson22pos  =  h/3  *  ssum22pos;  simpson23pos  =  h/3  *  ssum23pos; 

simpsonl2neg  =  h/3  *  ssuml2neg;  simpsonl3neg  =  h/3  *  ssuml3neg; 

simpson22neg  =  h/3  *  ssum22neg;  simpson23neg  =  h/3  *  ssum23neg; 

probmass2pos  =  l/pi*simpsonl2pos  -  l/pi*simpson22pos; 
probmass3pos  =  l/pi*simpsonl3pos  -  l/pi*simpson23pos; 

%  the  probability  mass  function  is  assumed  to  be  zero  for  the  final 
%  value  of  Z  and  thereafter;  probmass...  is  eq.  (29)  for  a  anu  b  values 
%  above . 

probmass2pos ( data+1 )  =  0;  probmass3pos (data+1 )  =  0; 
probmass2neg  =  l/pi*simpson!2neg  -  l/pi*simpson22neg; 
probmass3neg  =  l/pi*simpsonl3neg  -  l/pi*simpson23neg; 
probmass2neg(data+l )  =  0;  probmass3neg(data+l )  =  0; 

%  this  shows  how  close  the  sum  of  the  pmfs  is  to  1. 

disp( 'summation  error  (1/2)  pos  =  ');  disp( l-sum(probmass2pos ) ) ; 

disp( 'summation  error  (1/3)  pos  =  ');  disp( l-sum(probmass3pos ) ) ; 

disp( 'summation  error  (1/2)  neg  =  ');  disp( l-sum(probmass2neg) ) ; 

disp( 'summation  error  (1/3)  neg  =  ');  disp( l-sum(probmass3neg) ) ; 

stairs (d2 ,probmass2pos ) ;  pause  %  plot  the  pmf  in  stairs  mode, 
stairs ( d2 , probmass2neg) ;  pause 
stairs (d3 , probmass 3pos ) ;  pause 
stairs (d3 , probmass 3neg) ; 

D2(g)  =  sum ( sqrt( probmass 2pos  .*  probmass2neg) ) ;  %  this  is  the 

%  D-parameter 

D3(g)  =  sum ( sqrt( probmass 3pos  .*  probmass 3neg) ) ;  %  eq.  (30). 

%  the  following  formulas  are  approximations  to  the  probability  error 
%  bounds  for  rate-1/2  and  rate-1/3  convolutional  codes  —  eqs .  (31,32). 

bound2  (g)  =  36*D2(g)/'10  +  211*D2(g)~12  +  1404*D2  (g)  "14  ... 

+  11633*D2(g) "16  +  77433*02 (g) "18  +  502690*D2 (g) "20; 

bound 3 ( g)  =  ll*D3(g)"18  +  32*D3(g)"20  +  195*D3(g)~22  +  564*D3(g)"24  ... 
+  1473*D3(g) "26; 

disp( 'error  bound  for  rate-1/2  code  =  ');  disp(bound2 (g) ) ; 
disp( 'error  bound  for  rate- 1/3  code  =  ');  disp(bound3 (g) ) ; 
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